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A wide class of physical distributions appears to follow the (/-Gaussian form, which plays the 
role of attractor according to a g-generalized Central Limit Theorem, where a g-generalized Fourier 
transform plays an important role. We introduce here a method which determines a distribution 
from the knowledge of its q- Fourier transform and some supplementary information. This procedure 
involves a recently g-generalized representation of the Dirac delta and the class of functions on which 
it acts. The present method conveniently extends the inverse of the standard Fourier transform, 
and is therefore expected to be very useful in the study of many complex systems. 

I. INTRODUCTION 

Nonextensive statistical mechanics [H, Hj], a current generalization of the Boltzmann-Gibbs (BG) theory, is actively 
studied in diverse areas of physics and other sciences. This theory is based on a nonadditive (though extensive Q) 
entropy characterized by the index q, such that q = 1 recovers the standard BG entropy. It has been applied in 
systems such as cold atoms in dissipative optical lattices Q, dusty plasmas Q, trapped ions spin-glasses [7|, 
turbulence in the heliosheath j8|, self-organized criticality [9( , high-energy experiments at LHC/CMS/CERN Qjl and 
RHIC/PHENIX/Brookhaven low-dimensional dissipative maps [l2[, finance (ljjj . among others. 

This new theory can be advantageously based on ^-generalizations of standard mathematical concepts, such as 
the logarithm and exponential functions, addition and multiplication, Fourier transform (FT) and the Central Limit 
Theorem (CLT) H. Recently, plane waves, and the representation of the Dirac delta in plane waves have been 
generalized as well Some of these generalizations open the door to interesting aspects. For instance, a generic 

analytical expression for the inverse q-FT for arbitrary functions and any value of q does not exist (lij . The focus on 
this fact and related questions should be relevant for various applications in physics (e.g., field theory and condensed 
matter physics), engineering (e.g., image and signal processing), and mathematics where the standard FT and its 
inverse play a crucial role. 

We show in this paper that, in the 1 < q < 2 particular case and for non- negative functions (e.g., probability 
distributions), it is possible, by using special information that we shall detail later on, to obtain a biunivocal relation 
between the function and its q-FT. This property, not to be confused with an inverse q-FT, is not relevant to the 
proof in [lij of the g-generalized CLT, which yields g-Gaussian attractors (defined here below). It should be also clear 
that such an inverse g-FT is by no means necessary for the existence of attractors, which can be proved [20j without 
recourse to this (generically nonlinear) integral transform. 

The g-FT of a non- negative integrable function f(x) of the real variable x is defined as [l4j 

/•OO 

F q [f}(0= f(x)ef^^"- 1 dx (l<g<3), (1) 



where e z q denotes the principal value of [1 + (1 — q)z\ 1 ^ 1 ^ q ^ > (ef = e z ). Furthermore, the g-Fourier transform can also 
be defined for q < 1 (see |2l|). 

Another important connection concerns the Dirac delta. The 5 q distribution is defined as fl5j 

= 5— :/ (!<<?< 2), (2) 

For g = 1, this distribution is the usual plane- wave representation of Dirac delta. Also, the above integral corresponds, 
for an arbitrary g, to the g-FT of the constant function f(x) = 1. Wide families of functions exist |15l - fl8| such that 
the 5 q distribution behaves like the Dirac delta for g ^ 1 (see also section liTIj) . 



II. A METHOD ENABLING THE INVERSION OF THE g-FOURIER TRANSFORM 

Let 1 < g < 2 and f(x) be a non-negative piecewise continuous function of the real variable x, whose support will 
be denoted by supp/. Then, for each y S supp/, we define the function f^ v \x) = f(x + y) of the real variable x 



2 



(notice that f^°\x) = f(x)). The g-Fourier transform of f^ v \x) is given by 

/OO /"OO 
f(v)( x)e **V«H*)]<- 1 dx= / f( x + y )e^f^ q - 1 d x . (3) 
-OO J — OO 

Using the change of variables z = x + we have that 

J —oo -/supp / 

hence 

/• OO /»00 /• 

^[/ (y) ](C,y)^= / / fMep'-vM'^dzdt. (4) 

) J — oo J supp / 

Assuming that the function f{x) is such that it is allowed to commute the integral operators in the RHS of this 
equation, we have that 

F q \fW]{£,y)dZ= f f(z) r e «(»-»)I/wr 1 de < fa = -^- / f(z)S q ((z-y)[f(z)Y- 1 ) dz. 

' ./supp/ J — oo Q supp / 

Assuming also that the function f(x) belongs to the class of functions for which the 5 q distribution behaves like the 
Dirac delta. Then we have 



F 9 \f<» ) K&v)d£ = 7 ^ [ f(z)5((z-y){f(z)}i- 1 ) dz . (5) 

'supp / 



2-q 



Let us consider first the case in which supp / is a finite union of disjoint closed intervals, since in this case supp / 
has a boundary. Then supp / = Ufci -^i where Ig — [ag, be] with ag < be, and the Ig are mutually disjoint. Then 



where we used the property of the Dirac delta given in Appendix IA1 since fix) > for any x G supp /. As y G supp / 
has been fixed, there exists only one £q such that y G [ae , be ]- Then 



Therefore, 



f(v)=\^l F q [fM](t,y)dty , (6) 



77T 



DC 



where 7 = 2 for any interior point y 6 supp /, and 7 = 1 for any boundary point y G supp /. 

Let us now consider the case in which supp / = R. We straightforwardly obtain from Eq. ([5]) that 

F q [f^Kt v) r m§^r = ^imr q , 

2-qJ-oo [f{y)] q 2- 9 

which yields in the (7 = 2)-particular case of Eq. ((6]). 

If supp / is bounded either from above or from below, then Eq. (|6]) can be obtained from Eq. ([5]) following a 
procedure similar to the one described on the above lines. 

a. Particular case. Let us remark that, in the q — > 1 limit, we have that 

-1 />oo -1 y* 00 />oo 

f(y) = lT FlfW](Z,y)<%=— / f{x + y)e* x dxd(i. 
27r J -00 27r J -00 J -00 

Using the change of variables z = x + y, we have that 

-1 poo />oo -1 />oo />oo -1 />oo 

f(y) = — / f(z)e^ z -yUzdZ=- e-*»/ f(z)e*'dzd£=— e*" F[f](£) <% > (?) 
2n J -00 J -00 ^ J-00 J -00 ^ J-00 
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FIG. 1: Representation of G3/2,i(a;). The continuous line corresponds to the analytical expression of the function; the dots 
were obtained by handling numerically Eqs. (fjrj and JSJ. 

which is the well-known expression of the inverse Fourier transform. 

If we attempt to follow for q ^ 1 the above lines, we immediately find the difficulty that, in general, e q +b ^ e q e q . 
Therefore, no generic explicit expression analogous to (JJJ has been found for the inverse q-Fourier transform. 

Eq. © says that we can obtain the function / if we know the g-Fourier transforms of f{x) and of all its translations 
in the x-axis, which we illustrate next. 

Example 1. Let q < 3, j3 > and iel. A g-Gaussian distribution is given (see for instance [1, [3) by 



where C q is a normalization constant, in the sense that C q = c q x dx. Let us now consider 63/2,1(2;), then 



(8) 



G 



3/2,1 



(y) 



1 Z" 00 



(9) 



for 1 < q < 2. The case q = 3/2 was handled numerically (see Fig. [T]), since the g-Fourier transform of an arbitrary 
translation of G 3 / 2} i(x) could not be obtained analitically, even Eq. (J6]) being an analytical result. We also verified 
that Eq. (|9]) remains true for some other values of q. 



Example 2. Let q e (1, 2), A > 0, x € R and 

fA,g{x) = 



q-2 

\x\ 1- 1 - A 



\x\ 1- 1 



2(q-l) 

l + ( g _l)(|a;|l^_^ ' 2 



(10) 



when A < | a; |(9- 2 )/(9- 1 ) , and zero otherwise (see 19]). This function is constructed such that fo. q (x) — G q .i(x) and 
its q-FT does not depend on A. Thus, given only F q [fA, q ](0 = Fq[G q ,i](£) we cannot determine the original function 
fA,q- Nevertheless, Eq. (j6|) states that, if we know the q- Fourier transform of fA,q(x) and of all its translations in the 
x-axis, then we can determine fA,q(x) (see Fig. [2]). 



III. 5 a DISTRIBUTION AS THE DIRAC DELTA 



Let us now discuss a family of functions for which the 8 q distribution behaves like the Dirac delta. 
Theorem 1. Let q 6 [1, 2), and f(x) be a function of the real variable x which satisfies the following conditions: 

W JZo\f( x )\ dx <°°- 



4 



0.6 



0.4 



0.2 - 



01 




-1.0 



-0.5 



0.5 



1.0 



FIG. 2: Representation of fi t 5/i(x). The continuous line corresponds to the analytical expression of the function; the dots were 
obtained by handling numerically Eqs. (fTJ) and ([6]). For all values of x G ( — 1, 1) we have used 7 = 2 in Eq. JBJ, whereas for 
x = ±1 we have used 7=1. 



(") ri™(fc)|dfc<oo. 
Then, the 8 q distribution acts as the Dirac delta, i.e., 



f(x)S q {x)dx = f(0). 



Proof. If q = 1, we have nothing to add to what is already available in the literature. If q € (1,2), we can use the 
Gamma representation of the g-exponential, namely (l6| . 



1 



9-1 



e-H& e *(9-l)C«* dt . 



(11) 



hence 



Notice that 
1 



OO pCO 



f{x)S x dx 



f{x)e-H^e^ q ~ 1 ^ xt 



OO /»oo 



(A) 



f{x)e-Hyie l{q - 1 ^ xt dtdx. 



\f(x)\dx J e"H^dt< I \f{x)\ dx < 00, 



where we have used condition (i). Then, from Fubini's theorem, we can permute the integral operators, hence 



/OO -1 pOD pOO -f p 
f(x)ef dx = — ^ / e-Hf^ / f(xy^-^ xt dx dt = — ^ / 



r (^) 

We must now prove that 



e-H^F[f]{{q-m)dt. 



OO />OG 



2?r 



00 J — OO 



f(xYl*dxdt=- /(0) 

2-q 



(12) 



First, we should note that using the change of variables £ = (q — l)£i, we have that 



2-q 



OO /■ OO /■ OO 

/ e-*t53|jr[/](( ? _ i)^)| dlidt = / e~*t- 
J-00 Jo 

1 



m/KOI 



dt 



q-iJo 



e'H 1 ^ dt 



(q - l)t 
|F[/](C)|dC<oo, 
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where we used condition (ii). Then, using again Fubini's theorem as well as £ = (q — l)£t, we have that 



/oo poo -i poo pc 

/ f(x)ef x dxd£ = -. ^ / e-H^dt / 

-ooJ-oo (g- l)r ( ^j-J JO — ' 

2T/(0)r(fE?) _2tt/(0) 



^[/](C)dC = 



2tt/(0) 



- i)r (^t) 7 ° 



e'H^ dt 



(9-D r (^r) 



2-g ' 



which completes the proof. 



□ 



As a corollary, it follows from theorem [T] that the S q distribution behaves like the Dirac delta for the family of 
functions D ai b t \(x) = {ax 2 + b) x of the real variable x if a, b > and A < —1/2. Indeed, we can straightforwardly 
verify that 



7T/2 

|A, AA (a;)|da:oc / (cos^r 2 ^. 

-tt/2 



< OO . 



On the other hand (see [22] (3.771.2)) 



F[D aAX ](k) = 2 1+ \l — 




r(-A) 



where K n (x) is the modified Bessel function of the second kind. We can verify, using the asymptotic expansion of 
K n (x) (see (22] (8.451.6)), that F[Z? 0) 6 i >](fc) is a positive integrable function. As a particular case, it follows that the 
8q distribution with 1 < q < 2 behaves like the Dirac delta for g'-Gaussians with 1 < q' < 3. 

We analyzed numerically the behaviour of the S q distribution for the functions D a ^,\{x) with A > 0. Let us remark 
that these functions behave smoothly at x = 0, and diverge like a power law for |a;| —> oo. It results that there exists 
<Zmax(A) € [1, 2) such that the S q distribution behaves like the Dirac delta for any 1 < q < g ma x(A) (see Fig. [3]). 
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FIG. 3: For 1 < q < q , ma x(A), the 8 q distribution behaves like the Dirac delta for functions of the type -D a ,6,A(^), even when 
A > 0. The dots have been calculated numerically. The inset suggests that g max ~ 1 + l/A (A — > oo). 



IV. CONCLUSIONS 



Summarizing, we have shown that it is possible to determine the original distribution from the knowledge of the 
associated g-FT's for all possible translations y of its argument, in contrast with the knowledge of only the y = 
case, as shown in 19]. The present results clarify some questions that emerge in connection with a large number 
of experimental, observational and computational results in physical systems which do suggest the emergence of 
q-Gaussian and q-exponential distributions (0413 and others). 

The exact mathematical conditions under which the exchange performed in Eq. is legitimate remains as an 
interesting open point. Our examples reinforce however that it should be admissible for most physical cases. 
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Appendix A: A property of the Dirac delta 

If P(x) is a diffcrcntiable function of the real variable x, and P(x) has a finite number of roots, namely x±, X2, ■ ■ ■ , x n , 
such that P'(xj) ^ 0. Then 

in the space of continuous functions of the real variable x. 
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